Nonlocal Cooper pair Splitting in a pSn Junction 
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Perfect Cooper pair splitting is proposed, based on crossed Andreev reflection (CAR) in a p-type 
semiconductor-superconductor-n-type semiconductor (pSn) junction. The ideal splitting is caused by the en- 
ergy filtering that is enforced by the bandstructure of the electrodes. The pSn junction is modeled by the 
Bogoliubov-de Gennes equations and an extension of the Blonder-Tinkham-Klapwijk theory beyond the An- 
dreev approximation. Despite a large momentum mismatch, the CAR current is predicted to be large. The 
proposed straightforward experimental design and the 100% degree of pureness of the nonlocal current open the 
way to pSn- structures as high quality sources of entanglement. 



Spatially separated entangled electron pairs arise in hybrid 
normal-metal-superconductor structures. Andreev reflection 
(AR) is the conversion of an electron into a hole at the in- 
terface between a normal-metal or semiconductor and a su- 
perconductor |l|]. In nonlocal (or crossed) Andreev reflection 
(CAR) the conversion is over two electrodes maintaining the 
singlet state J^S]. This makes CAR a promising source of 
locally separated entangled electrons and building block for 
solid state Bell inequality experiments, quantum computation 
and quantum teleportation ||4j,|5|]. 

The Bell inequality can only be violated when the CAR 
fraction is larger than l/y/2 of the total current Despite 
intensive investigation, the CAR fraction is usually small due 
to competing processes yiW- Aside from local AR, tunneling 
of a particle from one lead to another can occur. Elastic cotun- 
neling (EC) is to lowest order in tunneling amplitude equal in 
magnitude and opposite in sign to CAR, resulting in a vanish- 
ing nonlocal conductance [9]. Including higher order terms, 
which become important in more transparent junctions, un- 
fortunately provides EC to be the dominant process iflOll - 

Several proposals have been put forward to enhance the 
CAR current. Using ferromagnetic-halfmetals (F) as leads 
can result in dominant CAR in an antiparallel magnetiza- 
tion alignment [3], though spin entanglement is then ques- 
tionable. Pure nonlocal Cooper pair splitting is predicted in 
a superconductor-topological insulator structure 111 111 , but the 
fabrication will be challenging. Both the electromagnetic en- 
vironment 111 211 and a change in the density of states (DOS) of 
a superconductor due to an ac bias II 1311 can result in dominant 
CAR, though the influences are expected to be small. A larger 
effect is predicted by Cayssol M1411 in an n-type graphene- 
superconductor-/?-type graphene junction. However, full can- 
cellation of EC and AR is only at precise biasing to the Dirac 
point and at a small range of the energy spectrum so that the 
current is not completely carried by CAR. Optimization of 
CAR may also be realized by using the Coulomb interaction 
0]. Recent experiments indicated great potential of using the 
energy to discriminate CAR ||8|], although the splitting effi- 
ciency is yet small. 

Here, we propose a strategy for ideal 100% nonlocal 
Cooper pair splitting, with no contributions from AR and EC, 
in a relatively straightforward device. Making use of the en- 



ergy difference between the incoming electron and the An- 
dreev reflected hole, in combination with bandstructure im- 
posed forbidden energies, enables the cancellation of both 
AR and EC individually, while still having a significant CAR 
probability. This idea is shown in Fig. [TJ The asymmetry of 
the bandstructure allows only CAR to occur since particles 
due to EC and AR will end up in forbidden states in the band 
structure, the band alignment being tunable by proper gating. 

A wide range of materials are suitable as electrodes. Ex- 
amples are nanowires, where the bandgap can be tuned by 
the length of the wire, minigap semiconductors such as bi- 
layer graphene, and narrow band semiconductors in general. 
Impurity bands can be used when there is significant density 
of states at an energy ±A from the semiconducting bandgap, 
generally in the case of low doping concentrations. Energy 
bands or levels that arise from quantum confinement in gen- 




FIG. 1 . (a) NSN structure with gates attached to the semiconductor. 
By applying gate-voltages, the semiconductor bands in each elec- 
trode can be tuned with respect to the Fermi energy. A bias voltage 
applied at the left NS interface results in a nonlocal conductance at 
the second interface, (b) In standard NSN structures an incident elec- 
tron can result in normal reflection (R), AR, EC, and CAR. (c) In an 
nSp junction it is possible to have perfect nonlocal Cooper pair split- 
ting. 
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eral can be used, as long as gapped energy regions exist that 
prohibit AR and EC. The proposed type of energy filtering is 
not only of use as Andreev entangler, but can also serve as 
energy beam splitter in, for example, FSF devices. This opens 
up an alternative route towards Bell inequality experiments; 
the spin may be employed to split the Cooper pairs, since the 
energy can be utilized for the read out. 

We model the NSN system by extending the classical 
Blonder- Tinkham-Klapwijk model [15] to three dimensions 
and by including the second interface. The BTK model has 
been used in modeling SNS and FSF structures where the di- 
mension of the sandwiched layer is close to £ ifl^l . In the 
present case of a pSn junction, the need to model in three di- 
mensions stems from the large Fermi momentum mismatch 
between a semiconductor and a superconductor. Because of 
the conservation of momentum parallel to the interfaces, a 
critical angle between momentum and interface normal ex- 
ists above which no transfer can take place. Each of the AR, 
CAR and EC probabilities is characterized by an energy and 
angle dependent effective barrier. Since excitation energies 
are comparable to or larger than the Fermi energy of the semi- 
conductor in our system, we will go beyond the Andreev ap- 
proximation that takes all momenta equal. 

We describe the pSn structure shown in Fig.Q]with the time 
independent Bogoliubov-de Gennes equations given by 



for "J in the regions N±, S, and N2, then becomes 



H(r) 
A*(r) 



A(r) 
-H(r) 



*(r) = EV(r). 



(1) 



^(r) = (u, v) T is the wavefunction in Nambu (electron-hole) 
space. We assume that A(r) = in the normal regions N\ 
(z < 0) and N2 (z > d), and A(r) = Ao in the superconduct- 
ing region S (0 < z < d), d being the superconductor width. 
The use of these rigid boundary conditions is warranted by the 
large Fermi momentum mismatch across the interfaces, which 
effectively reduces the coupling between the layers. A spec- 
ular barrier is included at the interfaces, resulting in U (r) = 
HiS(z) + H2S(d — z). Inside a region we assume an isotropic 
bandstructure. Incorporating the conserved momentum com- 
ponent parallel to the interfaces into the Hamiltonian allows us 
to simplify Eq. (Q]) to a ID system with an effective ID Hamil- 



tonian, given by H(z) 
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■U(z) — /j,*, where 



9* 2m*(z) 9z 

H* = (A* - E p ) cos 2 6± T \/E 2 - A 2 sin 2 B±. The angle 0± 
is the angle between the direction of the electron (+) or hole 
(-) and the normal of the interface and can be found through 
the Snell-Descartes law sin 6t = rj~ sin where and 9 t are 
the respective incidence and transmission angle and r& is the 
ratio of the incoming and transmitted moment. For large r*fc, 
8 t ~ and particles in the superconductor travel normal to the 
interface, fi is the chemical potential and E p is the potential 
energy in each layer tuned by the gate voltages. Our ansatz 
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= 1 - 
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and 



= |(1 + ^ E2 E A ' 2 ) [15]. Using conservation 



of fell we find the moments in the z-direction in each layer 
given by fc ± = 



cos ( 



'± 



2m • 



E p ± VE 2 - A 2 ) 



. This 

implies a critical angle given by 9c{E) = arcsin(r^ 1 ). In a 
pSn junction with E p 7^ in the electrodes, CAR is enhanced 
since CAR has the lowest critical angle as compared to the 
other scatter processes lfl4ll . 

The system can be solved by applying the boundary 
conditions ty Nl 2 (0~ , d+ ) = * s (0 + ,d _ ) together with 

^^(o + ,<n- 4 ^ (o-/)=ig^.iM. 

S JV 1>2 

extended BTK model reproduces the results in NSN structures 
found previously by other models. The nonlocal conductance 
Gnl vanishes due to the cancellation of CAR by EC in the 
tunnel limit IM], while EC is dominant in transparent regimes 
ifioll . and the electrode separation distance dependence of Gnl 
is exponential. Charge imbalance is not taken into account in 
this model, but the CAR enhancement effects as described in 
this Letter also occur in the range E < A at temperatures 
T <C T c , where this effect is absent 11711 . 

In order to enhance the CAR current to 100%, p and n type 
semiconductors will be implemented now as the two elec- 
trodes. We investigate the generic example of a semiconduc- 
tor possessing a band gap Eb ^S> A, one valence band and 
one conduction band, each with parabolic dispersion. Even 
scattering processes at energies within the forbidden semicon- 
ducting band gap have nonzero probabilities at the interfaces, 
influencing other processes, and thus need to be taken into ac- 
count when solving the wave equations. Still, the wave func- 
tion of these particles decay exponentially in the electrodes 
and do not contribute to the current by themselves, since we 
consider electrode lengths much larger than £. 

Figure |2] shows the numerically obtained probabilities for 
perpendicular incidence, zero bias, transparent interfaces (the 
BTK barrier strength Z = = IU5IP and in a regime of 
large momentum mismatch. The probabilities are given by the 
absolute squared values of the prefactors in Ansatz Eq. (0. 
Normal reflection (R), EC, AR, and CAR probabilities are 
considered at the respective interfaces, while the quasiparticle 
(QP) states in the superconductor are taken far from the inter- 
face thereby vanishing at energies below A. In the regime of 
large momentum mismatch, the angle dependence below the 
critical angle and the bias and gating dependence of the prob- 



Po(E)\ 



E-eV-E c 



I cos 2 ((9), where 



abilities follow P{E) — , uliJ) 
Pq refers to the unbiased probability at perpendicular inci- 
dence. As an example, we consider Al as the superconductor, 
resulting in a large effective barrier originating from a large 
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FIG. 2. Energy dependence of the averaged reflection and transmis- 
sion probabilities. Z — and E c \ =fx= E V 2 = 5.8 x 10 4 A, result- 
ing in an effective energy dependent barrier. Consequently, the prob- 
ability (CAR) is significant and almost equal to (EC) for E<A. 
Because of the small superconducting width d = rr£, (AR) does not 
reach unity and (R) does not vanish at E—A. The inset shows the 
dependence of the probabilities on the width of the superconductor 
at E — i A; Fabry-Perot resonances occur on the Fermi wavelength 
scale. 



FIG. 3. Local (ii) and nonlocal (I2) current dependence on the 
bias voltage (Vi) for a fixed gate voltage eV g i = 0.6A (see upper 
inset). E cl + 0.6A =fj, = E v2 = 5.8 x 10 4 A, d = and Z = 0. 
The momentum mismatch results in an effective barrier. 7?o is the 
Sharvin resistance at eV = A. Negative biasing (eV<0) leads to a 
nonlocal current due to EC, whereas positive biasing (eV>0) results 
in CAR. AR is possible up to 0.6A, so that in the range 0.6A<eV / <A 
perfect Cooper pair splitting occurs. QP current appears at eV> A 
shown in the lower inset. 



ratio [i/A = 5.8 x 10 4 . Despite the effective barrier, CAR is 
found to have a considerable magnitude for d being close to 
the Bardeen-Cooper-Schrieffer coherence length £ = VF ^ . 
Averaging the probabilities by (P) = f^u^Xpg P{z)dz 
with Xf.s the Fermi wavelength in the superconductor is nec- 
essary in order to let Fabry-Perot resonances vanish, in prac- 
tice washed out by roughness. 

The current density in the electrodes is obtained from 

J d = Y E J 9 (k,±)e z /(k,±). 

k,(T,± 

Here, d is the dimension and Va the volume of the system, the 
sum ± is over electrons and holes, / is the nonequilibrium 
distribution and the charge current is defined by J g , e (ft) — 
f|lm [u*(v*)Vu(v)] . From this it follows that the CAR and 
EC current are opposite in sign, since the respective group 
velocities in N 2 , given in the ansatz Eq. (O, are opposite 
in sign. Finally, the current in an electrode becomes / = 
2^ / dEf J Q f dO cos 9 sin 6 Y,± |k|mJ|, where A is the 
cross-sectional area. The integration is over all energy modes 
that contribute to the tunneling, limited by the lowest DOS of 
the initial and final state for a certain process. Even nonideal 
pSn junctions with nonvanishing DOS in the bandgap or im- 
proper Fermi level aligning enhance CAR, since the DOS will 
be lowest for the AR and EC processes. The distribution func- 
tions are given by /1 = fo(E — eVjvi) — fo(E — eVs) and 
k,e{h) = fa{E - eVm) - f (E - emax[V s ; V N2 (-V N2 )]), 
with Jo (E) the Fermi distribution function. Positive or nega- 
tive biasing at the second electrode decreases the EC or CAR 
processes, respectively. Current flowing to the superconduc- 
tor is defined positive, so AR and CAR are positive in sign 
and EC negative. 



Figure [3] shows the IV characteristics for a pSn junction 
with fixed gate voltages, so that E c \ + 0.6A = // = E v2 - 
At negative bias across the first interface, the first electrode 
has available states above and below /i whereas the second 
electrode has only states below //. AR and EC are therefore 
possible, while CAR is prohibited. For positive bias voltages, 
direct electron transfer is no longer possible and the nonlocal 
current is carried by CAR only. AR is significantly reduced by 
the critical angle and limited DOS and totally vanishes above 
eV>0.6A, where available states for AR are absent. The de- 
vice works as a perfect Cooper pair splitter for 0.6A<eV< A. 
Above the gap a quasiparticle current appears in the supercon- 
ductor resulting in a lower splitting efficiency. 

Maximizing the bias regime in which CAR = 100% can be 
achieved by tuning the gate voltage such that always E c \ = 
fi = E v2 , irrespective of the bias voltage. This bias situation 
is shown in Fig. [4] For all positive bias voltages below A, 
CAR = 100%. AR and EC are forbidden due to the bandgap 
in the electrodes. The DOS for incoming electrons is equal 
to the DOS for outgoing holes and there is no critical angle 
lowering the CAR probability, since Ep = 0. Consequently, 
the nonlocal current is maximized to a relatively large value, 
being typically a few uA. 

We now address the question what materials would be best 
suited. The discrimination between AR, EC, and CAR lead- 
ing to 100% crossed Andreev reflection in a pSn junction is 
caused by: the forbidden band gap (which is our main effect 
leading to 100% CAR), the variation in density of states and 
the critical angle. The elasticity mandatory for these three 
effects is typically a less stringent condition than the super- 
conductor width being comparable to £, necessary to have 
a significant nonlocal current. With current nanolithogra- 
phy methods and using Al as superconductor these require- 
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FIG. 4. Local (7i) and nonlocal (J2) current dependence on the bias 
voltage (V2) for ideal gating, E c i =(M = E V 2 = 5.8 x 10 4 A regard- 
less of the bias voltage. d= |£ and Z = 0. The momentum mismatch 
results in an effective barrier. Ro is the Sharvin resistance at eV — A. 
For eV<0 no process is possible, but in the range 0<eV<A only 
CAR is possible and we observe pure entangled current I\ = l2- QP 
current appears at eV>A, lowering the CAR fraction. 



ments are easily fullfilled. For the semiconductors, InAs 
two-dimensional electron gases are ideal candidates, since no 
Schottky barrier is formed in contact to Al II 1811 . The effect of 
the critical angle may vanish in the diffusive limit, but perfect 
Cooper pair splitting due to the forbidden bandgap remains 
robust against disorder. Nb/InAs structures are therefore also 
suited. Even though £jv6 < Cai, the larger superconducting 
gap of Nb simplifies the band alignment and increases the 
magnitude of the nonlocal current. Schottky barriers reduce 
the nonlocal current, but the splitting is found to remain ideal 
when using a nonzero barrier strength. Al/GaAs heterostruc- 
tures are, therefore, suitable as well fl9tl . Finally, we mention 
that electronic gate controllable InAs nanowires have been 
contacted to Al with high interface transparency [20], mak- 
ing it an ideal system for entanglement experiments where a 
reduced number of propagating modes are required. 

In conclusion, we have proposed a pSn junction that can be 
used to prepare a pure Bell state by forward biasing and can 
act as a perfect Cooper pair splitter by reversed biasing, while 
having significant currents. 

Discussions with A. A. Golubov, B. C. Kaas, C. J. M. Ver- 
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